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ABSTRACT: Using gauge/gravity duality, central ultrarelativistic nucleus-nucleus col- 
lisions are modelled as collisions of shock waves in five-dimensional asymptotic AdS 
space. For early times after the collision, it is possible to analytically match the 
metric from the past to the future light-cone. This allows extraction of the pre- 
equilibrium energy-momentum tensor of the strongly coupled, large N gauge theory. 
For central collisions, this allows qualitative statements concerning the build-up of 
radial flow in AA and pA collisions. We find that the early-time radial flow buildup 
is identical to that expected from ideal hydrodynamics with an entropy density pro- 
portional to the square root of the product of the matter densities in the individual 
"nuclei". 
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1 Introduction 

A long-standing problem in the field of high energy nuclear collisions has been trying 
to understand the precise mechanism and timing of equilibration after the collision 
of the incoming nuclei. The success of the heavy-ion program at the Relativis- 
ts Heavy Ion Collider (RHIC) [1-4] and the Large Hadron Collider (LHC) [5-7] 
has further provided motivation to study this difficult regime, where one needs to 
describe the real time evolution of non-perturbative coupled quantum field theory. 
Previous work on the pre-equilibrium regime in nuclear collisions can broadly be clas- 
sified in two branches: weak-coupling (perturbative) techniques to solve Quantum- 
Chromodynamics (QCD) [8-16] and gauge/gravity duality to solve strongly coupled 
gauge theories different from QCD[17-29] . In the former setup, partice-like degrees 
of freedom are weakly coupled to the classical field background (such as that from the 
Color-Glass-Condensate) and are known to give rise to a kind of plasma instability 
[30] . At the time of writing, the actual equilibration (transition to hydrodynamic be- 
havior) has not been observed in realistic (e.g. longitudinally expanding) simulations, 
but lots of progress has been made in toy model systems. 

This work will follow the second, strong-coupling approach, where one permits 
oneself to trade the non-abelian gauge theory of physical interest (QCD) with a 
different (non-abelian) gauge theory N — 4 SYM. In this case, using the dictionary 
from gauge/gravity duality [31], one can map the energy-distribution of a fast moving 
nucleus to a gravitational shock wave in asymptotically AdS space in five-dimensions 
(cf. [20]). Using machinery from relativity, one can draw on previous work that has 
established the form of the line element for given (energy-) density profiles of nuclei, a 
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lower bound on the total entropy produced in the collision of two of the shock waves. 
Moreover, using techniques from numerical relativity on the collision of black holes 
in various dimensions, the actual equilibration of the gauge theory from a far-from 
equilibrium state to hydrodynamics has been observed in numerical simulations with 
a high degree of symmetry (the nuclei where assumed to be translationally invariant 
in the plane transverse to the collision axis) [25, 26]. However, these simulations 
could not provide information on the pre-equilibrium dynamics in the transverse 
plane, which is of interest because it could potentially lead to observable effects in 
nuclear collision experiments at RHIC and the LHC The aim of the present article is 
to provide the foundation to lift this shortcoming by providing the metric shortly after 
the head-on collision of two gravitational shock waves, including the extraction of the 
early-time gauge-theory energy-momentum tensor. Presumably, using these results 
in conjunction with recent numerical advances in solving the Einstein equations in 
AdS^ (cf. [32]) will allow observation of equilibration of the system including full 
transverse dynamics. 

2 Setup: Heavy-ion collisions as gravitational shock waves 

Boosting a charge to very high velocities, its energy-density distribution becomes 
highly singular. This is very similar to the case of boosting a mass to very high 
velocities. However, in the latter case it has been understood by Aichelburg and 
Sexl in the 1970's that a reasonable description of the energy-momentum tensor can 
be given by rescaling the mass with the boost factor [33] . By complete analogy, the 
energy-momentum tensor of a boosted charge p can be calculated analytically by 
means of a rescaling of the coupling constant [34] , arriving at a form of 



where light-cone coordinates x ± = ^| have been introduced and x± = (x, y) are the 
coordinates in the plane transverse to the boosting direction. 

Within gauge/gravity duality, it is known how to construct a strongly coupled 
M = 4 SYM configuration that has precisely this energy momentum tensor. The 
line element is given by 



where z is the coordinate parametrizing the fifth dimension in a space that is asymp- 
totically Anti-de-Sitter and 



+ + oc p(x ± )5(x + ) , 



ds 2 



2dx + dx~ + dx\ + dz 2 + dx +2 $(z, x_ L )5(x + ) 
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where k is the simulated number of degrees of freedom. Note that while for M = 4 

TV 2 

SYM k = it is easy to make the number of degrees of freedom more QCD-like 
by resetting k by hand, as done e.g. in [26]. 

This line element corresponds to a gravitational shock wave in AdS^ with a 
transverse profile that is governed by the density distribution p{x±). It is an exact 
solution to Einstein equations if [35] 



d 2 z -- z d z + dl 



®(z,x ± ) = 0, (2.2) 

which can be solved in Fourier-space to give [36] 

${z,k ± ) = c 2 {k ± )z 2 I 2 {zk ± ) , 

where I 2 is a modified Bessel function. Note that the second solution to (2.2) does 
not describe a M = 4 SYM field theory in Minkowski space and therefore is not 
allowed here. The function c 2 (k±) is related to the Fourier-transform of the profile 
function p(kj_) by 

ip(k ± ) 



c 2 (k±) 



One should note that one could supplement the gauge/gravity setup by intro- 
ducing sources (currents) in the bulk, in order to regulate the behavior of the profile 
function $(z, x±) for large values of z. However, it turns out that for the purpose of 
shock wave collisions, the results including these sources is identical to those without 
sources [37]. Therefore, we chose to not modify the original gauge/gravity setup and 
work with profile functions that may not have well defined z — > oo limits. Since in 
the process of the shock wave collision, a horizon will form at a finite value of z, the 
large z behavior is no longer relevant for the subsequent evolution anyways. 

As a model for a boosted nucleus, the form of p can be taken to be a Wood-Saxon 
distribution, a distribution calculated from the Color- Glass- Condensate model or any 
other favorable model. The present technique for the pre-equilibrium dynamics is 
not limited to any specific nuclear physics model for a boosted nucleus. 

For example, one can take $(z, k±) oc z 2 I 2 (zk±)K 2 (zok±), which in configuration 
space can be shown to be 

\x ± \) <xzq~ 3 2^1(3,5/2, 3, -1/q), 

where q = X±+ ^ Z J°^ and 2 F\ is a hypergeometric function that takes a particular 
simple form [20, 36]. This particular choice corresponds to a profile function p(x±) oc 

(4 + ^ 2 r 3 - 
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3 Methodology 



The first important point to address is generalizing the matching conditions for the 
collision of two shock waves with transverse profile p = const in Ref. [22] to arbitrary 
profiles p(x±). To perform the matching, it is advisable to transform the line ele- 
ment to a different iable form that does not contain a 5-function (Rosen coordinates). 
For a general transverse profile p(x±), this can be easily done by generalizing the 
corresponding case for spaces that are asymptotically Minkowski [38]. The relevant 
coordinate transformations for a single shock wave are 

1 u8 2 (u) • 1 

x + =u, x~ =v+-$6(u) + —±-!-d i $6 t3 d j $, x l = x 1 + -uO^S^d^ , (3.1) 

2 8 2 

where the coordinates x l = (x, y, z) were introduced. Superposing the metrics from 
shock wave one and shock wave two one obtains a line element that is valid before 
the collision: 

-2dudv + d&dxi (s kl H^H^ + 5 kl H^Hf - 8 i3 ) 

_ V { J _ / /g <j\ 

pre [z + \{u9{u)d z m + vO{v)dM 2) )] 2 

where u < 0, v < and 

Central Collisions While the following program seems straightforward for any 
(physically allowed) pair of source functions <$>@\ one can expect the calculation 
to be rather tedious. Therefore, in this article we limit ourselves to considering 
head-on collisions of azimuthally symmetric nuclei, for which azimuthal symmetry is 
unbroken in the future light cone. Hence one introduces new coordinates r, 

x = r cos 0, , y = r sin , z = z 

for which one finds in particular 1 

because $ = $(r, z). For better readability, it is convenient to introduce the notation 

The pre-collision line element ds pie then becomes 

-2dwcfo + dx l dx j Mf- e 
dst - 3 



r " [z + \ (u9(u)d z $W + v0{v)d z &° 1 1 ' 



(2) 



1 Note that the arguments of $ in Eq. (2.1) can be take as x l because the ^-function multiplying 
$ erases the distinction between these. Clearly, the same argument will not hold true if one had 
had "smeared" delta functions instead. 
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with 



Mr = i + ««(«)•&; + ((*^) 2 + (*;;;) 2 ) 

ac = i + + ^ + (*;;',) 2 ) + [« « «, i « 2] 



JMT = + !L ^-4>( 1 ;; (*m + *?i>) + [« « «, 1 ** 2] , (3.3) 

and all others vanishing. 

Matching The matching at the collision point is performed by making an ansatz 
for the line element such as 

ds 2 = ds 2 pTe + 9(u)9(v)ds 2 post (3.4) 

with ds 2 ie given in Eq. (3.2). The number of independent metric functions appearing 
in ds 2 ost can be found be with the following argument. First note that if the line 
element ds 2 was in Feffer man- Graham form, then this would fix the metric function 
g zz to be 1/z 2 and g uz , g vz , g rz , g^ z to vanish. The no n- vanishing components of the 
metric would be in a sub-matrix spanned by the coordinates u, v, r, <ft, so in principle 
there would be 10 independent components. Limiting myself to systems that have 
azimuthal symmetry implies that g u </>, gv</>, gr</> are vanishing, and hence one only has 
7 independent metric functions. For symmetric collisions (e.g. $i = $2); one may 
thus pose 



ds 2 ost = —2dudvM uv + y-dv 2 H — du 



uu 

V U 

.2 71 r 1 j j.2 



M„ 



+2dr(du + dv)M ur + dr 2 M„ + d<p 2 M H + dz 2 M zz + 2drdzM rz . (3.5) 

Asymmetric central collisions (to model for instance proton-nucleus (pA) collisions) 
will be treated separately below. Expanding around u = 0,v = 0, we demand 
that the Einstein equations are fulfilled across the the light-cone. In essence, this is 
equivalent to the procedure performed for the so-called Color Glass Condensate [39] , 
where one deals with Yang-Mills equations instead of Einstein equations and one 
matches gauge field configurations instead of the metric field. 

Close to the light-cone u ~ 0, v ~ 0, so one can make an ansatz such as 

M rr (u,v,r,z) = ( u + v )f^( r ,z)+uvf 2 ^r,z) + (u 2 + v 2 )f 22 (r,z) 

+uv(u + v)f%(r, z) + (u 3 + v 3 )f™(r, z) + uV/ r 4 r °(r, z) 
+uv(u 2 + v 2 )ft 2 (r } z) + (u A + v A )f^(r, *) + ..., (3.6) 
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and equivalently for the other metric functions M a ^. Note that absence of negative 
powers of u, v simply follows from the observation that these would give rise to highly 
singular terms such as 5'(u)/u in the Einstein Equations. 

In order to solve Einstein's equations on the light-cone (u — or v — 0), the 
coefficient for all singular functions has to vanish. Specifically, one finds that the 
condition that there are no S(u)/u, 5 2 {u) or 5'{u) terms in the u, v component of 
Einstein's equations immediately leads to 

./// oj 22 (.../- 0,# 0... 

meaning that the line element has to be continuous across the light cone u — 0, v — 0. 
Using this information, it is vastly more convenient to again change coordinates using 
the so-called Milne form 

, 1 M 

r = V2uv , f = - In - , 
2 v 

where 

—dr 2 g TT + T 2 d£ 2 g££ + dr 2 g rr + r^dcf^g^ + 2drdzg rz + 2rdrdrg Tr + dz 2 g zz 



ds< 



and metric functions 



g TT = 1 + K(r,£,r,z) , 

9tf = 1 + L(r,£,r,z) , 

g rr = Mr r e (l + H(T,£,r,z)) , 

9^ = M; r ;(l + F(T,Z,r,z)) , 

g rz = M? r z e + G(r, £, r, , 

g Tr = J ( T) ^ r , , 

^ 2 = Mf; e (l + M(r,e,r,^)). 



(3.7) 



(3.8) 



For early times r « 1, the seven metric coefficient functions J, K, L, H, F,G, M 
may be expanded in a Taylor series around r = and the coefficients of this taylor 
series are determined by solving the Einstein Equations order by order in r. The 
resulting line element ds 2 may then be brought into a more convenient form, such as 
Feffer man- Graham coordinates or Eddington-Finkelstein coordinates by a suitable 
coordinate transformation. 



4 Solution for central AA collisions 

Using the methodology outlined in the last section one can find a solution the case of a 
head-on collision of two shock waves with identical profile functions, $(i) = $( 2 ) = $. 
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This could be interpreted as a model for nucleus-nucleus (AA) collisions, such as Pb- 
Pb at LHC energies. In this case, the metric functions J, K, L, H, F, G, M can be 
expanded as 

K = T 2 k 2 o(r,z) + r 3 cosh £ k 31 (r, z) + r 4 fc 40 (r, z) + r 4 cosh 2£ fc 42 (r, 2;) + . . . , 

which corresponds to the expansion in Eq. (3.6) plus the additional knowledge that 
K must vanish for u = or v = (absence of singularities on the lightcone). 
Solving the Einstein Equations order by order in r one finds for example 

k 20 2 ($ 0,1 ) 2 $0'! ($!.0 - r$ 2 ' ) ($ 1,0 ) 2 $ 1 >°$ 2 >° ($i.i) 2 + ($2,0)2 

^20 — — H 



20 



3 3z 2 2rz 6r 2 6r 6 

($0,1) 2 ( $ 1,1) 2 + ($ 2,0)2 



Az 2 



/20 "~^ 2_ + ^ 2 ~' (4 ' 1} 
3$o,i $i,o $ i,i 

920 = ^2z 2T~ ' 

17 ($0,1)2 3^0,1^1,0 ($1,0)2 ($1,1) 2 3$0,1$2,0 $1,0 $ 2,0 ($2,0)2 

4z 2 rz 2r 2 2 z r 2 ' 

while the expression for j 2 o is too lengthy to be reproduced here. Note that the 
coefficient function K is only determined as a constraint at higher orders. Specifically, 
one find that its early-time, near-boundary expansion is given by 

With the post-collision metric known, one would like to extract information 
about the boundary energy-momentum tensor T^ v where \i = 0, 1, 2, 3 and x M = 
(t, x, This is most easily achieved by rewriting the line element in Fefferman- 
Graham form , 

1 9 u-v dx^ doc ~\~ dz z Tnydx^doc ^ z hu,v doc^doc , * 

ds = + ^- + > , (4.2) 

z z z A z £ 

n=0 

where the boundary metric is assumed to be flat: g^ v = diag(— 1, 1, l,r 2 ). To bring 
the line element Eq. (3.7) into the Fefferman-Graham form, one uses the coordinate 
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transformation 

oo 

T = T FG + ^2 tn(TFG, O^G, r FG )z p + ^ n » 

?1=0 

OO 

£ = £FG + ^ e »( r re, ^FGi r F G)z A F G 2n , 
n=0 

oo 

r = r FG + ^2 s n(r FG , £ FG , r FG )z F + Q n , 

n=0 

oo 

2 = ^FG + ^2 a n( T FG, &G, r FG )z F + G 2n , 



n=0 

^FG, 



(4.3) 



with coefficients t n , e n , s n , a n which are determined order by order in z. One finds 
that to lowest order the required coordinate transformations are 

COsh^G$°' 4 4 

t0 = 247^ + ° ( } ' 

smh^ ! + 

S °" 2471 ( } ' 

a o = - rreCO ^ G$O ' 4 + 0(r 5 ). (4.4) 

Since at z — z FG = we have t> g = r, the resulting boundary energy-momentum 
tensor may be written in terms of the boundary coordinates x M . Solving Einstein's 
equations up to (including) 0(r 3 ) one finds 

288 v ' 

K = 96~ + ( } ' 

($0,4\2 / $0,4^1,4 \ 4 

T " = TiT 7 * + ( 6 " ^T- " 3 * M * 2 ' 4 ) TT52 + 0(r5) ' 

T H = + (2 (*") 2 - 3*-^- - 4>».^) + C(r=) , 

$0,4^1,4 3 

288 V ; ' 

T ?r = G(r 5 ) , 

T €r = 0(r 5 ) . (4.5) 
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Up to the order given, this energy-momentum tensor is traceless and covariantly 
conserved, V^T 7 "' = 0. In fluid dynamics, it is custom to decompose the energy- 
momentum tensor using the energy density e and the fluid velocity with = 
— 1. Even for non-equilibrium (dissipative) fluids, these are completely defined in 
terms of the eigenvalues and eigenvectors of the energy momentum tensor 

UfjT^ u = -eu u . 

Ignoring for the moment the fact that the solution Eq. (4.5) is not that of a fluid, 
one may still ask what values of e and the energy-momentum tensor corresponds 
to if one pretended it was that of a fluid. One finds 

f(T>0> 4 ^ 2 <J>1,4 

'=^88^' - r = -^ + <^ 3 )> ^ = 0(r). (4.6) 

This finding implies that there is a radial flow u r building up that is proportional to 
the gradient of the transverse charge distribution of the nucleus, 

d T u' = - drH 3 Pir)) , (4.7) 

which should be compared to the evolution expected for ideal hydrodynamics with 
conformal equation of state (c 2 . = 1/3) [41]: 

<9Xdeal hydro = ~$Pr h S , (4.8) 

where s is the initial entropy density. Clearly, the evolution found for here for the 
early-time pre-equilibrium radial flow is identical to that from linear ideal hydro- 
dynamics, provided one identifies the square root of the initial overlap distribution 
a/ p 2 (r) with the entropy density s(r). The same is not true for the extracted energy 
density e, which according to linear hydrodynamics should decrease rather than in- 
crease. In any case, the present calculation provides a concrete example for a far-from 
equilibrium system evolution with a flow profile identical to that expected from ideal 
hydrodynamics. If a similar phenomenon were to happen for anisotropic flow, this 
would have important implications for the attempt to use experimental anisotropic 
flow measurements to extract the viscosity coefficient of hot QCD matter. 



5 Solution for central pA collisions 

The above strategy to find the line element after the collision of two shock waves 
may be generalized to the case of asymmetric collisions ($1 ^ $2), which may be 
taken to represent a model for the proton-nucleus collisions (pA). In this case, there 
is no longer a symmetry u -<->■ v and hence it is adviseable to slightly generalize the 
ansatz for the post-collision line element ds pos t (3.5). Matching predicts that the line 
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element is again continuous across the light-cone, so one may switch again to Milne 
coordinates r, £. We find that the line element 



ds 2 



x [-dr 2 g TT + T 2 d^,(g^ + 2g^ r dr + 1g iz dz) 
+dr(g rr dr + 2g rz dz + 2rg rr dr) + r 2 d(f) 2 g^ + dz 2 g zz ] , (5.1] 



is suitable for obtaining a solution to the Einstein Equations, even though (or maybe 
because) it contains two redundant metric functions. Following the strategy outlined 
in section 4, one can solve the Einstein Equations order by order in a power series ex- 
pansion in t. The resulting solution may then be transformed to Fefferman-Graham 
coordinates and one finds the following result for the energy-momentum tensor: 



N 0,4^0,4 



288 V ) 

$°' 4 $°> 4 



_2 Jl 

96 

=,0,4^0,4 



T„ = ^^r 2 + 0(r 4 ) , 

$0,4^0,4 
rp ^1 ^2 J2J2 | 

T H = 7-7-, r r + 0(r 



144 

1,4 
2 

576 



T rr = ^ 2 1 2 -^— + C(r 5 ) 



TV = 0(r 5 ) , 

$1.4 $ 0,4 _ ^0,4^1,4 

TV = 1 2 1Q2 1 2 r 4 + G(r 5 ) . (5.2) 

As it should, this result corresponds to the form for AA collisions when $x = $ 2 - 
Note that in this case one obtains a pre-equilibrium radial flow 

_ ($?< 4 ^' 4 + ^' 4 $°' 4 ) r 

uT ~ 6$r$F ' (5 - 3) 

which again can be interpreted as ideal hydrodynamic flow buildup (4.8) for an 
entropy density 

s(r) oc \/pi(r)p 2 (r) , 

with pi,2(^) the transverse charge density distribution of nucleus 1,2, respectively. 



6 Summary and Conclusions 

In this article, we studied the central collision of gravitational shock waves in asymp- 
totic AdS$ spacetimes. For early times after the collision, we find an explicit result 
for the post-collision line element. The novel aspect about our study is that we 
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allowed the shock waves to have arbitrary profiles with azimuthal symmetry in the 
plane transverse to the collision axis, thus generalizing the result of Ref. [22]. Via 
gauge/gravity duality, we are able to interpret our result as the early-time energy 
momentum tensor generated by the collision of two "nuclei" in A^ = 4 SYM. 
Our three most important findings are: 

1. It is possible to generalize the methods developed in cf. [22] to the case of in- 
cident shock waves with less symmetry, thereby indicating that the technique 
should also be applicable to gravitational shock waves with no special symme- 
tries in the transverse plane. 

2. The resulting early-time energy-momentum tensor reflects the fact that the 
system is initially far from equilibrium. For instance, the effective longitudinal 
pressure is negative. Our result for the pre-equilibrium energy-momentum 
tensor implies matter flow in the (transverse) radial direction. This serves as a 
concrete example for the generation of far-from equilibrium flow in high energy 
"nuclear" collisions. 

3. The build-up of this radial flow is identical to that expected from ideal hy- 
drodynamics with an entropy density proportional to the square root of the 
product of the charge density of the individual shock waves. This last result 
is somewhat unexpected because the system, being far from equilibrium, does 
not evolve according to ideal hydrodynamics as a whole; only the radial flow 
buildup does. 

We also studied the implications for non-symmetric collisions as a model of 
proton-nucleus (pA) collisions. In this case, all three of the above points also apply. 
In particular, our result implies strong (as compared to AA collisions) early-time 
radial flow as a consequence of the proton charge density falling of much more steeply 
in the radial direction than that for a heavy nucleus. This could possibly explain 
some of the experimental findings in proton-lead and proton-gold collisions at high 
energies. 

The biggest limitation of the present study is that it is only applicable at very 
early times where a power series expansion converges. However, there is a good 
chance that these early time results may be used as input for a subsequent numerical 
solution of the Einstein Equations (e.g. from Ref. [32]) to obtain full results for all 
subsequent times, as in Ref. [26]. We intend to pursue this direction in a future 
publication. 

As another application of the present work, we want to emphasize that there 
seems to be no obstacle in generalizing the present results to the case of shock wave 
collisions with arbitrary transverse profiles p(x±), such as provided by the PHOBOS 
Monte-Carlo Glauber model [40], which results in a lumpy distribution for p. In this 
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case, subsequent results for the post-collision dynamics would have to be averaged 
over many configurations (events), yet from the experience gained in the present 
work, also this case seems to be feasible using known techniques. 
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